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Abstract 

We develop a bosonization approach for one-dimensional models based on 
Bethe ansatz equations. The operator formalism of the exact soluble models 
in the low energy limit provides a systematic method to calculate the asymp- 
totic correlation functions. As examples with and without internal degrees 
of freedom, the Calogero-Sutherland (C-S) model and the repulsive Hubbard 
model are considered respectively. We verify that the low energy behavior 
of the C-S model is controlled by two classes of c = 1 conformal field theo- 
ries, depending on whether the C-S interactions are among bosons or among 
fermions. For the Hubbard model, we show the explicit charge-spin separation 
at low energy for arbitrary U > 0. The low energy behavior of the system is 
described by the (semi-) direct product of two independent Virasoro algebras 
with c = 1. 
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I. INTRODUCTION 



Recently it has attracted much attention to understand strongly correlated fermion sys- 
tems in low dimensions {d < 2). Thanks to the special behaviors of one- dimensional (1-d) 
systems, we have better knowledge for them than those in higher dimensions, with the 
help from some non-perturbative methods. There is a number of 1-d systems called ex- 
actly soluble, whose spectra are exactly given by the Bethe ansatz equations. However, it 
is often difficult to calculate the correlation functions by using the Bethe ansatz. On the 
other hand, the bosonization approach provides a useful tool to compute the asymptotic 
correlation functions [0] in certain cases. Therefore, it is important to understand the inter- 
relationship between exact solutions and bosonization. According to Haldane the Bethe 
ansatz soluble models can be characterized by 1-d Luttinger liquids. Moreover, after renor- 
malization, only forward scatterings are important in the Luttinger liquids. These imply 
that there is a harmonic fluid description |^ to exact soluble models in the low energy limit, 
i.e., it should be possible to bosonize these models, at least in the low energy regime. How- 
ever, how to do bosonization, especially for the case of the systems with internal degrees of 
freedom, has never been shown explicitly before in the literature. In this paper we will show 
how to bosonize the exact soluble models at low energy by starting from the Bethe ansatz.. 
The bonus of doing so is a short cut for deriving the low temperature thermodynamics and 
static correlation functions. We will take two nontrivial models, the Calogero-Sutherland 
(C-S) model and the Hubbard model as examples. 

The C-S model not only presents a beautiful example for soluble 1-d many-body 
model with long range interactions, but is also closely related to physical phenomena such as 
the edge excitations in the fractional quantum Hall effect 0. It also provides a realization 
of some basic physical concept such as exclusion statistics MM. The correlation functions 



of the C-S model can be exactly calculated in a wide range of the coupling constant |10,11 



Also, Ha [|T^ proposed a bosonized anyon effective theory at low energy, and found that 
the asymptotic correlation functions of the effective theory agree with those of exact results. 
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This raises a question of whether there is an exphcit way to derive the bosonization and 
show this agreement. In our recent work [jl2|, a bosonization approach was developed for 
1-d gas with exclusion statistics, which can be modified to deal with this problem. We 
will apply this approach to the bosonization of the C-S model, as a warm-up for treating 
models with internal degrees of freedom such as the Hubbard model. A by-product of this 
exercise is the proof that at low energy the bosonic and fermionic C-S models give rise to two 
different c = 1 conformal field theories (CFT's), because of different selection rules for the 
quantum numbers. To our knowledge, this result is new in the literature. We also present 
the correlation functions respectively in each case. 

The Hubbard model is known to have many remarkable properties, and has attracted 
a lot of attentions since the two dimensional Hubbard model was believed to be the right 
model for the high Tc-superconductivity about ten years ago. Although the 1-d Hubbard 
model we studied here may not directly relate to two-dimensional physics, there are still 
good reasons to investigate it for gaining insights into high Tc-superconductivity. After 
Lieb-Wu's Bethe ansatz solution of the 1-d Hubbard model [0, a lot of efforts have been 
spent in this model. We list here some relevant developments. As is well-known, there may 
be string solutions with complex rapidities in some Bethe ansatz soluble models, which may 
affect the low energy physics of the system [l^. Fortunately, for the Hubbard model, it has 



been shown that the complex rapidities offer only the states in the upper Hubbard band and, 
therefore, do not contribute to the low energy physics [|1^. So we are allowed to consider 
only the real rapidities in our treatment. Spin-charge separation is a fascinating property of 
the Hubbard model. It has been shown that in the strong-coupling limit, the ground-state 



wave function given by Lieb-Wu's solution is indeed of the spin-charge separation form [jT6 
It was also checked from the exclusion statistics point of view that the low-lying states are 
also of the spin-charge separation form for large U Ren and Anderson assumed 

the separation for arbitrary coupling in the low-energy limit and proposed a bosonized 
effective theory of the model. The critical exponents of the effective theory agreed with the 
numerical calculations |1T6| , |T8| and the phase shift argument fl^. The correlation functions 



of the Hubbard model were also discussed by finite-size scaling considerations [|T9[. It has 
been proven that there are two independent c = 1 Virasoro algebras which describe the 
critical theory of the Hubbard model. It was also asserted that the low energy behaviors of 
the system at less than half filling could be characterized by a (semi-) direct product of two 
Virasoro algebras |]TU|. 

In this paper, we consider the repulsive Hubbard model at less than half filling. The low 
energy effective theory of the Hubbard model is not Lorentz invariant in 1 + 1 dimensions, 
since there are two different Fermi vectors. This generally leads to the gap-less charge- 
and spin-density waves having different velocities. In this sense, the spin-charge separation 
is expected to appear not only at large U but also at low energy at arbitrary coupling, 
as has been assumed in ||T^. Starting form the Bethe ansatz solutions, we will explicitly 
demonstrate this property in the low temperature thermodynamic limit. We show that the 
free energy of the model at the low temperature is of the form F{T)/ L = F{0)/L — — 
which implies that the critical behavior of the system can be characterized by the 
direct product of two independent c = 1 Virasoro algebras and coincides with the results 
of finite-size scaling at T = Meanwhile, by generalizing the bosonization approach 
(without internal degrees of freedom) to the Hubbard model which incorporates spin degrees 
of freedom, we arrive at an effective theory by bosonizing the thermodynamic limit of the 
original theory, resulting in the bosonized effective theory proposed by Ren and Anderson. 
Then the calculations of the asymptotic correlation functions become systematic. 

This paper is organized as follows: In the next section, we present bosonization of the 
C-S model. In Sec. HI, some relevant results about the Hubbard model are reviewed. Then 
in Sec. IV, the thermodynamic potential of the Hubbard model (and then free energy) at low 
temperature is derived from the Bethe Ansatz equations and the spin-charge separation at 
low energy is demonstrated. In Sec. V the bosonization of the Hubbard model is given, and 
the single-particle correlation functions are systematically calculated for various excitations. 
The last section is devoted to the conclusions. 



4 



II. CALOGERO-SUTHERLAND MODEL 



The C-S model we consider here is described by the A^-body Hamiltonian 

N Q2 

H=-T.^ + T.U{x.-x,), (2.1) 
i=i '-'^i i<j 

with 

U{x) =g + ^L)-' = g^, sm-\^), (2.2) 

where L is the size of 1-d ring. As L ^ oo, the interaction potential ( p.2|) U{x) —>■ g/x^. 
The ground state is given by 

i^B^ = n I sin 7r(xi - x,)/L|\ (2.3) 

or 

^F^ = X{{z,-z,)'Xlz,'^-'^'\ (2.4) 

i<j k 

depending on the statistics of the particles, bosons or fermions. Here z = exp(i27ra;/L) and 
A = [(1 + 2gY^'^ + l]/2. It is obvious that for A = 1(5' = 0) the former is the ground state of 
hard-core bosons and the latter that of free fermions. 

Moreover, the model can be exactly solved by using the asymptotic Bethe ansatz 
or the Jack Polynomials [1^. The periodic boundary conditions give rise to the spectrum 



determined by the following equations for pseudomomenta kf. 



Lh = 2Tih + 7r(l - A) ^sgn(A;, - h), (2.5) 



where 2 1 



{(A^ + l)/2 mod(l), for fermion case, 
(2.6) 
integer, for boson case 

These selection rules have been discussed by Kawakami and Yang, and for bosons and 



fermions the behavior of the momentum distribution is quite different ||2^ . We will see that 



it is the difference in the selection rules that gives rise to two different classes of c = 1 CFT 
for low energy effective field theory. The total energy E and the total momentum P are, 
respectively, 

E=j:ki p=J:k^=J:Mk). (2.7) 

i i i 

where Pi{k) = 27r/j/L 

The thermodynamics of the theory has been fully discussed by Sutherland p. Here 
we would like to examine its relations to CFT, Luttinger liquid and bosonization. Some 
relevant topics have been partially considered by several authors before The 



thermodynamics of the C-S model can also be formulated in terms of an ideal excluson gas 
(lEG) [1,0. The thermodynamic potential is known to be given by 



n _ T 
L ~ 



/oo 
dkln{l + w{k,T)-^), (2.8) 
-co 



with the function w{k,T) = pa{k)/p{k), where Pa{k) and p{k) are the hole and particle 
densities, satisfying 

w{k, T)^[l + w{k, T)]'-^ = e(^'-'^)/^. (2.9) 

The particle density is determined by 

p{k, T)(l + w{k, T)) = - (A - l)p{k, T). (2.10) 

Ztx 

In the ground state, there is a (pseudo-) Fermi surface kp, such that p{k) = l/2n\ for 
|A;| < kp and p{k) = for \k\ > kp. Then the Fermi momentum is given by kp = nXdo, and 
the ground state energy and momentum by Eq/L = tt^A^c/q/S, Pq = 0, where Jq = Nq/L. 

Now let us examine possible excitations in the model. First there are density fluctuations 
due to particle-hole excitations, i.e., the sound waves with velocity Vs = vp = 2kp (see 
below). Moreover, by adding extra M particles to the ground state one can create particle 
excitations, and by Galileo boost a persistent current. It is easy to verify that the velocities 
of these three classes of elementary excitations in the model satisfy a fundamental relation, 
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Vs = y/VNVj, that Haldane years ago used to characterize the Luttinger hquid |2^. Indeed, 
shifting Nq to N = Nq + M, the change in the ground state energy is 6iEq = 'K{\kp)M'^ , 
while a persistent current, created by the boost of the Fermi sea k —>■ k + nJ/No, leads to the 
energy shift 62EQ = 7T{kp/X)J'^. Therefore the total changes in energy and in momentum, 
due to charge and current excitations, are 

6Eo={7r/2L){vNM^ + vjJ^), 

6Po = TT{do + M/L)J, (2.11) 

with 

Vn = VsX, Vj = Vs/X, Vs = yJVNVj. (2-12) 



These coincide with the well-known relations |^ in the Luttinger liquid theory, if we identify 
A with Haldane's controlling parameter exp(— 2(y9). It is easy to check that the selection rule 
^!6D can be rewritten as [pO|] 

{Mmod(2), for fermion case, 
(2.13) 
even, for boson case 



It has been claimed from finite-size scaling pO[ that the low energy behavior of the C- 
S model is controlled by c = 1 CFT. A review of this point will help to understand the 
behavior of the Hubbard model. The central charge c of a conformal invariant system is 



related to thermodynamic observable by finite-size scaling Generally, the free energy 
of the system is given by 

L L ~ 6vs' ^ ' 

For the C-S model, we consider the low energy expansion of the thermodynamic potential 
( p.8|) . Following Yang and Yang p7|^, we introduce the dressed energy e{k,T) by writing 



ti;(A;,T) = e^('='^)/^. (2.15) 

And (^) reads 



e{k, T) = k^~fi-T{l- A) ln(l + e'^^^'^)/^). (2.16) 

Because there is no singularity in e{k,T) at T = 0, the zero temperature dressed energy is 
given by 

r - kl)/x, \k\ < kp, 

e{k) = {' (2.17) 
l> P - kl, \k\ > kp. 

At low energy, one can consider only the excitation around the Fermi surface and the 
thermodynamic potential is given by |^ 



n{T) 1 r^F ^^^^^^ _T d5ppHl + e-l^^'^^^)'^)!/^), (2.18) 



L 271 J-kp n J-s 

where 5 is a cut-ofF and the first term on the right hand side of the last equality is recognized 
as Q{0)/L. 6pF is the deformation of the Fermi surface, i.e., Spf/^tt = p{kF)6kF- It is known 
that e{k) = VgSpF near the Fermi point pf- After a little calculation, ( p.l^ ) reads 

fi(T) Q{0) _ nT' 

~L ^-"6^' ^^■^^> 

which implies that the theory is actually cut-off independent at low temperature. Notice 
that F = Q ~ pN. Only the particle-hole excitations near the Fermi surface contribute 
to thermal excitations, leading to N{T) — N{0) = 0; this can be checked by an explicit 
calculation in terms of the definition ( |2.1(j| ) of p{k,T). Thus, we have 

F{T) F(0) Q{T) Q{0) irT^ 



L L L L 6vs 

Comparing this to (|2.14]) , we see c = 1 for the C-S model. 



(2.20) 



The central charge of a conformal invariant system can also be given by the finite-size 
scaling in the spatial direction [^. In Appendix A, we confirm c = 1 for the C-S model 
from such a calculation, with some subtleties not noticed before in the literature. Now let us 
work out explicitly the effective theory at low energy with c = 1 conformal invariance. For 
this purpose, we will develop a bosonization approach. Our starting point is the observation 
that the grand partition function Zq, corresponding to the thermodynamic potential ( p. 81) , 
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is of the form of that for a system of ideal fermions with a comphcated, T-dependent energy 
dispersion given by the dressed energy: 

ZG = Y[{l + e-<''^y^). (2.21) 

k 

This fermion representation is not very useful, because of the implicit T-dependence in the 
dressed energy. Nevertheless, as we have seen in ref. in the low-T limit the dressed 

energy is effectively T-independent: Namely we have e(fc,T) = e{k) + 0(Te~''^*-'^-"/"^), as 
a result of ( |2.16| ). Hence, the grand partition function can be obtained from an effective 
Hamiltonian, given by 



H,s = J2<k)cick, (2.22) 

k 



where are fermion creation operators. 

Another simplification in the low-T limit is that we need to consider only low-energy 
excitations near the right and left Fermi points \k\ = ^kp, which are completely decoupled 
and have a linearized energy dispersion 

f ±VF(k^^kF), \k\ > kp, 
e±{k) = ^ ^ ^ I I (2.23) 

±vp{k^kp)/\, \k\ < kp- 
We note the 'refraction' at k = ±kp. In spite of this peculiarity, we have succeeded in 



bosonizing the effective Hamiltonian as follows ||T2[: The density fluctuation operator at 
k kp is constructed as follows 

k>kp k<kp—Xq 

+ E -t^,.,,/^^ (2-24) 

kp-\q<k<kp ^ 

for g > 0. A similar density operator p^~^ can also be defined at A; ~ —kp, 

k<—kp k>~kp+\q 

+ E ■■c\_^_o,: (2.25) 



—kp+\q>k>—kp 



— kp-q 



Within the Tomonaga approximation [p^], in which commutators are taken to be their 



ground-state expectation value, we obtain 

[P?\pfV^?V, [H^,P^t^] = ±v,qp^^\ (2.26) 
which describe 1-d free phonons with the sound velocity Vg = vp. Introducing normal- 



ized boson annihilation operators bq = y27r/gLp^+) and bq = y27r/gLp^ ^'1', the bosonized 
Hamiltonian satisfying ( |2.26| ) is given by 

He = vAY: q{h\bq + b\bq) + \j[\M' + (2.27) 

g>0 ZL \ 

The bosonized total momentum operator, corresponding to the fermionized P = 
J2kP{k) clck, is 

^ = E ^(blbq - b\bq) + 7r( Jo + M/L) J. (2.28) 

In the coordinate-space formulation, the normalized density field p(x) is given by p{x) = 
Pr{.x) +Pl(x): 



Pr{x) = ^ + E ^q/2nLX{e^''^bq + (2.29) 

^ q>0 



and Pl{x) is similarly constructed from bq and 6j, 

Pl{x) = ^ + E ^/q/2nLXie-'''% + e''^%), (2.30) 

where Mji ^ are given by M = M/j + M^. The boson field (/'(x), which is conjugated to p{x) 
and satisfies [0(x),p(x')] = i6{x — x'), is 0(a;) = 4>r{x) + 0l(x) with 

Mx) = ^ + ^ + ^Y: ^nX/2qLie^'^%q - e'^^-fet), 

(2.31) 

with J = Jr+Jl- Here M and J are operators with integer eigenvalues obeying the selection 
rule ( |2.13| ), and 0o is an angular variable conjugated to M: [0o,^] = i- The Hamiltonian 
( p.271 )becomes 
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Hb = ^ r dx [U{xf + {d,X{x)n (2.32) 

ZTT JO 

where U{x) = n\^/^p{x) and X{x) = \^^/^(p{x). With X{x,t) = e'"^X{x)e-'"\ the La- 
grangian density reads 

£ = — (9„X(x,t)9"X(x,t). (2.33) 
2tx 



We recognize in no ambiguity that C is the Lagrangian of a c = 1 CFT [0. Since 4>q 
is an angular variable, there is a hidden invariance in the theory under — > + 27r. The 
field X is thus said to be "compactified" on a circle, with a radius that is determined by the 
coupling constant: 

X ~ X + 2'kR, = 1 /A. (2.34) 

(Incidentally, A can be identified with the exclusion statistics parameter [|l^ of the quasi- 
particle excitations in the C-S model.) States V^[X]|0) or operators V[X] are allowed only 
if they respect this invariance, so the quantum numbers of quasiparticles are strongly con- 
strained. For the boson case in (|2.13|) , the zero- mode term in (|2.27| ) can be written as 
(7r/2L)[MV-R^ + 4:R'^D'^], where D = J/2 = Mmod(l) are integers. This means that the 
Hamiltonian ( p.27| ) in this case has duality 



R^1/2R, M ^ D, (2.35) 



which is just the famous duality in the usual c = 1 CFT ^D[. For the fermion case in ( ^.13p , 
the duality of ( p. 271) reads 



R ^ l/R, M ^ J, (2.36) 
which corresponds to the so called excluson particle-hole duality mentioned in 

A ^ 1/A. (2.37) 
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Moreover, the partition function in the fermion case (in the low-T hmit) can be rewritten 
as Z = Trfi\q^oq^o^^ where q = e™"^'^. The zeroth generators of the Virasoro algebra are 

Hr =VsY.^b% Hl =VsY.qh% (2.38) 

The constraint J = M mod (2) makes the spectrum and duahty relation different from those 
in the c = 1 CFT compactified on a circle. Indeed, the fermion case corresponds to the c = 1 
CFT compactified on an interval (or an orbifold) S^/Z2; see ||30| for more details. 



To show the difference between two CFT's, we check the primary fields in each theory. 
The primary fields of both CFT's are given by 

The charge-1 operators correspond to M = 1. Therefore, charge-1 operators are given by 
M = 1 and J = 2m with integer m for the boson case, while J = 2m + 1 for the fermion 
case. Hence, the charge-1 primary fields in the boson case are bosons, and those in the 
fermion case are fermions: 

oo 
m=— oo 

. gi(Ai/2_2m/AV2)x^(x+) . ^ ^2.40) 

and 



^^(x,t) =p(x)i/' ^ e'°™ : e 

m=— oo 

^i(Ai/2_(2m+l)/Ai/2)Xi(a;+) 



oo 

iOm . ^j(Ai/2 + (2m+l)/Al/2)Xfl(a;_) . 



. gi(AV--(2m+iVA^'-jAi(,a;+j . ^ (2.41) 

where the hermitian, constant- valued operators Om satisfy [Om-, Om'] = «7r(m — m'), which 
give rise to the Klein factor necessary for the correct commutation relations for \E'b(_f) and 

It is clear that the correlation functions of the bosons are different from those of the 
fermions, e.g., for m = 0, one has respectively 
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(*B(x)*yO)) ~ (2.42) 

and 

(^j.(a;)^U0)) ^ cos A;i.a;. (2.43) 

In summary, we have bosonized the C-S model and shown that the low energy behavior is 
controlled by c = 1 CFT. Depending on whether the C-S interactions refer to those between 
boson or between fermions, we have two different classes of c = 1 CFT, governed by different 
selection rules for quasi-particle quantum numbers. To our knowledge, the appearance of 
two different classes of c = 1 CFT's in the low energy limit of the C-S model was not notified 
before in the literature. 

In addition to the above charge- 1 operators, there are more allowed operators, the so- 
called anyons or exclusons, in the theory constructed as follows: 

^l^^{x) =: ^lj(a;)e*^'^'(^«(^)-^^(^)), (2.44) 

for boson case and 

*t,^^(x) =: ^t,(^)gi(AV2_i/Ai/2)(XH(x)-x,(.))^ (2.45) 
for fermion case. Those operators satisfy 

'^{{x)^l{x') - e''^^*»"(^-^')^^(a;')^^(a;) = 
for X /^x'. The multi-sector density operator for exclusons is 

p{x) = '^i{x)'^x{x) = *;j(^)(x)*s(F)(x) 

= p{x)^^ : exp{t2m[Xn{x) - Xl{x)]/X'^^} : (2.46) 

The correlation functions for density fluctuations and for anyons (or exclusons) are the 
same in the above-mentioned two classes of c = ICFT's: 



j2 



{p{x,t)p{0,0))^d, 



1 1 1' 

^ ^ (27rJo)'A 1^4 ^ 



^ ^/"' [x+X-l^VA cos(27rJoma;) 
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(2.47) 



G{x,t;X) = {¥^{x,t)^x{0,0)) 

-4 E ^^^^^^TX-^^'^'-''''^^'^^ (2-48) 

with Am and Bm regularization-dependent constants. We notice that the correlation func- 
tions (|2.48|) coincide with the asymptotic ones calculated in [|l^] for the C-S model. 



III. THE LIEB-WU SOLUTION AND RELEVANT DEVELOPMENTS FOR 

HUBBARD MODEL 

In the previous section, we achieved bosonization of the C-S model which has no internal 
degree of freedom. The bosonization of 1-d models with internal degrees of freedom such as 
spin is more complicated, because there are backward scatterings and umklapp scatterings. 
We will take the Hubbard model as an example to investigate. In this section, we first review 
Lieb-Wu's Bethe ansatz solution and other developments which are relevant in the coming 
sections. 

The general form of the Hubbard Hamiltonian reads 

{ij),a i i i 

where the summation {ij) runs over nearest-neighbors. We have set the electrons on a 1-d 
ring with size L. This model has been exactly solved by Lieb and Wu |Q, and the spectrum 
of the model is shown to be determined by the following Bethe ansatz equations, 

Lfc. = 2vr/. + y:2tan-^'(^^"^--^"\ 

E 2 tan- ^^^^-^^^^) = 2.J^ + 2 tan^ ^i^£L^, (3.2) 

where N (N^) is the number of electrons (with down-spin) and u = U/t. The total energy 
and momentum are related to the pseudomomenta ki by 

N N 

E = -2tJ2 cos ki + fiN + h{N^ ), 

i=i 2 



N 

(3.3) 



i=l ^ \ i a / 
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Although the variation of rapidity seems to be related to only the density variation 
of spin-down electrons in the A-space, it reflects the spinon density fluctuation in the real 
space, because of the antiferromagnetic ground state. Hence, the spin wave rapidity A^ 
characterizes the spinon dynamics of the system. The quantum numbers Jj and Ja obey the 
selection rules 

/. = ^mod(l), J, = + V od(l). (3.4) 

We would like to discuss the low energy behavior of the model with repulsive interaction 
{U > 0) at less than half filling. We have two kinds of gap-less excitations. Fortunately, it 
has been proven that the excitations with complex k and A have gaps of the order u = U/t 



15|. We can only consider the real k and A in the low energy limit. It has been shown that 



the finite-size corrections to the thermodynamic potential at T = are |19 



^l{0) - l^oo(O) = - — {v, + Vs) + 0(1/L2), (3.5) 

which implies that the low energy effective theory of the Hubbard model is not Lorentz 
invariant. There are two Fermi velocities Vc and Vg, which correspond to the charge- and 
spin- wave velocities respectively, and hence four Fermi points ±kp and ±Xf. The finite-size 
scaling ( p.5|) also shows that the low energy spectrum can be characterized by two c = 1 
Virasoro algebras. To derive ( p.5|) , it is important to realize that not only the thermodynamic 
potential (or the total energy) is scaled but also the distributions of the charge and spin- 
wave. In the next section, we will give ( |3.5|) another consistent check from finite-size scaling 
in the temperature direction. 

Again, we proceed to consider the thermodynamic limit. The zero temperature behavior 
has been well-understood recently ||T9|. Around the Fermi points, there are three kinds 
of excitations. The first kind is charge- and spin-wave fluctuations in real space, which 
corresponds to non-zero modes in k-X space. We shall leave those to the coming sections. 
Other two kinds are corresponding to zero modes in k-X space. One of them is created by 
adding extra number of particles to the ground state. The other is persistent currents. Both 
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of excitations shift the ground state energy. It is known that the ground state energy per 
site is given by 

1 [''F 

^0 = ^cik), (3.6) 

where Ec is the zero-temperature dressed energy obeying 

zvr J-Xp 
1 r^F 

SsiX) = 4 +7;- dk cos kKi{\ - sin k)ec{k) 

-1. d\'K2i\-\')ssi\'). (3.7) 

ZTT J-Xp 

The bare energies e^^^ = ~ h/2 — 2t cos k And ef^ = h. And the functions Ki and K2 are 
defined by 



The solutions of the integral equations (|3.7|) define the energy bands. So, we can determine 
the kp and A^? by 

£e(±A;F) = 0, es{±XF) = 0. (3.9) 
Alternatively, we can fix the values of kp and A^ through the particle numbers: 

nc = — = / dkpc{k), 



L J—kp 

Ni r^F 

ni = -r= dXp^iX), (3.10) 

J ~Xp 

where the distributions of charge- and spin-wave, p{k) and P|(A), are given by a set of 
integral equations: 

Pcik) = TT + ^ c/Airi(sinA; - A)p|(A), 
1 r^'F 



Pil-^) = / dkKi{\ - smk)pc{k) 

dX'K2{X-X')pi{X'). (3.11) 
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Generally speaking, the zero mode excitations can be thought as the fluctuations of the 
Fermi surface. The corresponding shift to the ground state energy is given by 



6E, = ^i5k,r + ^i6X,r, 



dkl 



(3.12) 



where the fact ( |3.9| ) insures there are no cross derivatives. It is easy to see d'^Eo/dhj 



and d^Eo/dXl 



Vs, where 



27rp|(Ai.)' 



(3.13) 



2ti pcikpY 

are the charge- and spin-wave velocities. 

Adding extra particles M and M| to the ground state leads to a deformation of the Fermi 
surface given by 



/ \ / dkp dkp \ 



6 m Xp 



dN dNi 
dXp dXp 



( m\ 



(3.14) 



V dN dNi I 

For the persistent currents, which correspond to the quantum number shifts 1^1 + Dj 
and J ^ J + Dj, we have the Fermi surface deformation as follows 



ojkp ' 



6jXf 



' dkp dkp \ / \ 



dDi dDj 
dXp dXp 



(3.15) 



\dDi dDj J 

In general, the matrices in ( |3.14| ) and (|3.15|) are related to the dressed charge matrix of 
the theory |T9[. For the case of zero magnetic field, h = 0, which we are interested in, we 
have 



^Mkp — 
5jkp 



2^ L '■ 
2L'' 



SmXp - 
5jXp 



1 

2^2 L 



where = M, Ms = - M^, = 2Di + Dj and Js = -Dj. And ^ 
the function ^{x) determined by integral equation 



(3.16) 
^(sin kp) with 
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1 rs'mkp 

^{x) = l + — / K{x~x')ax'), 

^TT J — sinkp 



p—\uj\u/A+iuix 

Kix) = I -duo. (3.17) 



In particular, we have 

p^{k) = ^ + / dk'K{sm k - sin k')p,{k'), (3.18) 

ZTT 2n J-kp 



e^{k) =e^^\k) + — dk' cosk' K{smk - smk')e^{k'), (3.19) 
27r J-kp 

Hence, the zero mode part of the excitation spectrum at /i = 0, according to our above 
discussion, is 

<5Eo = (7r/2L)(t;,,^M2 + v.^jJl) + {t: /2L){vs,nMI + ), (3.20) 
where the velocity relations are given by 

Vc{s),N = Vc{s)K{s), Vc{s),J = Vc{s)/K{s), (3.21) 

with Ac = and = 1/2. The energy shift ( |3.2q ) and the velocity relations ( p.21| ), for 
each component (spin or charge), resemble those in the C-S model, ( p.ll| ) and ( p.l2| ). Also, 
it is easy to show the momentum shift caused by zero modes is given by 

SPo = 7r(noe + MJL)J, + ttM^Js/L, (3.22) 

which can be thought of as the two-component generalization of the momentum shift in 

The selection rules (|3.4| ) imply that the excitation quantum numbers obey the following 



selection rules [19 



Di = ^ ^ mod(l), Dj = ymod(l). (3.23) 



For Jc and Jg, the fact = Msmod(l) and ( |3.23| ) imply that 



J, = ^mod(l), Js = ^mod(l). (3.24) 
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IV. LOW TEMPERATURE BEHAVIOR AND SPIN-CHARGE SEPARATION 



In this section, we would like to present the thermodynamics of the Hubbard model at 
low temperature and show spin-charge separation for arbitrary U > 0. Also, we will confirm 
the finite-size scaling result ( p.5|) . 



A. Statistics and thermodynamics 

Recently an exclusion statistics description has been developed for the Bethe ansatz 
soluble models P, p3| , [l4| . For the Hubbard model, the statistics matrix is given by 



g'^'^{k,k')=6{k-k'), 



TTM 1 + 16(sin k — Xy/u'^ ' 



gii{X,X')=6{X-X') + - ^ 



nl + A{X- X'y/u^' 



g^%X, k) = -1^^^—^^. (4.1) 
vr 1 + 16(A — sm ky /w^ 

The equation for g^-^ shows the mutual statistics between different spin-wave rapidities, while 
g'^^ and g^'^ give the mutual statistics between k and A. However, the constraint ( |3.9| ) means 
that the low energy excitations near the Fermi surface have no mutual statistics between 
states described by k and A. 

We showed before that for Bethe ansatz soluble models with no internal degree of free- 
dom, one can use an effective statistics parameter to characterizes the statistics of low energy 
excitations, which can be read off from the Luttinger-liquid velocity relations [|l2l. For the 



present case, the general velocity relation ( p. 211) seems to suggest that there is no effective 



mutual statistics and there are two effective statistics parameters, Ac and A<j, characterizing 
the statistics of low energy excitations. Below, we will see that indeed there are two types 
of anyonic (or exclusonic) excitations in the theory. 
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If only the real k and A are taken into account, the thermodynamic potential can be 
written down in terms of the principle of exclusion statistics and the statistics matrix 
([4.1[). It has been given in p^ : 

n T 

L 

where the dressed energy ec{k,T), which is the finite temperature generalization of ( |3.7|) , is 
given by 

J" POO 



-L rdfcln[l + e-^=(^'^)/^], (4.2) 



ec{k,T) = ef^ - — dXKi (sin A; - A) 
X ln[l + e-^^("'^)/^], 
( A, T) = r dk cos kKi{\- sin k) 



X ln[l + e 



ec{k,T)/Ti 



+ / dX'K2{\ - A') ln[l + e-^^(^''^)/^]. (4.3) 
27r J-oo 



The dressed energy es{X,T) does not explicitly appear in ( [4.2|) since the number of bare 
single particle states for spinon is zero. 

Corresponding to the zero temperature dressed energy ( p.l9| ) at h = 0, its finite temper- 
ature counterpart is 

e,{k,T)=e^^ 

-— r dk' COS k'K{sin k - sin k') ln[l + e~'^^^^^\ (4.4) 
27r J-TT 

Similarly, the charge- and spin-wave densities at finite temperature are also given by a 
set of integral equations 

p,(fc,T) = [l + e-('='^)/^]-^(;l 

L /TT 

+ 



/ d\Ki{smk-\)p^{\T) 



27T 

p,(A,T) = [l + e-(^'^)/^]-^ 



1 

— / dkKi{\ — smk)pc{k) 

27T J-TT 

dX'K,{X-X')p^{X') 



1 

Yn J- 



oo 



(4.5) 
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B. Finite-size scaling and spin-charge separation at low temperature 



In this subsection, we would like to do finite-size scaling at finite temperature, to show 
two new results. First, we want to confirm the finite-size scaling ( p.5|) at finite temperature, 
as promised above. On the other hand, the spin-charge separation in the Hubbard model 
has been previously shown only in the strong coupling limit for a given temperature |TB|,in 



Here we want to show the spin-charge separation at low temperature for arbitrary U > 
which was assumed in the Luttinger liquid interpretation of the Hubbard model []T7[. 

The thermodynamic potential ( [4.2|) looks like that of a single-component system ( p. 81) 



because the dressed energy es{X,T) does not explicitly appear. However, it is necessary to 
emphasize that edk^T) is coupled to Ss^X^T) through ( ^^ . Even at T = 0, the dressed 
energies are still coupled each other through ( |3.7|) . To see the finite-size scaling in T, we fix 
N = and = No/2. As shown in Appendix B, T) = e^^s){k{X),0) + O(^), 

i.e., one has 

e,{k, T) = e,{k) + e,{k, T) + OiT'/u^), 

BsiX, T) = es{X) + £s(A, T) + 0(rVz/2). (4.6) 

According to (|4.3[), ec{k,T) and es{X,T) can be determined by the following integral equa- 
tions 

e,{k,T) = --——K,ismk-XF)+ / —6,{X,T), 
e,{X,T) = -£^K,{X-XF) 

— -i^'ifsin kp — X) cos kp 

6e'^{kF) 

— COS kKi{sm k — X)ec{k, T) 
-kp 2tt 

r^t dX' 

- —K,{X-X')e,{X,T). (4.7) 
Solving the integral equations, one finds that 
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where 



So, we have 



/e = 27rp,(M - 1, A = 27rp,(A^). (4.9) 



2n J-kp ' 6ff. V 2-11 pikp)) 6fo 



-fcF 6fc \ 2TTp{kp)/ UVs 

Because we have fixed the particle numbers, the free energy can be expanded as 
T r , T r dk 



2tx 



-/iiVo-^(iV^-iVj 

+— / dk€,{k,T), (4.11) 



where F(0)/L = so ((§]§)). Using ( ^TTUQ , we have 



F(T) F(0) vrT^ vrT^ 



(4.12) 



L L 6vc Qvg 

In fact, this result is valid for arbitrary N = Nq + M and A^^ = A^o/2 + M| if M < Nq and 
< Aro/2. 

Finally, we get the free energy in the low temperature limit, consistent with (|3.5| ) and 
resulting in two c = 1 Virasoro algebras. Combining ( 4.12| ) and ( 3.20 ), we see that at 



h = 0, the contributions of excitations, including both zero modes and non-zero modes, to 
the thermodynamic potential are separated into spin and charge parts. In this sense, we can 
say that there is a spin-charge separation in the Hubbard model at low temperature. 
In the low-T limit, the thermodynamic potential can be rewritten as 

L L TX 

2T 

I,{kp,T)27rp,{kF), (4.13) 

71 
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where 



Is{Xf,T)= I " (iAln(l + e-^=W/^), 

rkp+S 

IcikF,T)= dk\n{l + e"'^^''^/^. (4.14) 

Jkp 

The physical momenta corresponding to the charge and spin excitations are defined by 

dp = 2TTpc{k)dk, dp = 2TTpc{X)dX. (4-15) 
Finally, the thermodynamic potential in the low-T limit can be expressed as 

-e,(k,±)/T)^ 



- r rf5p,ln(l + e-l"^('^P^)l/^), (4.16) 

TV J-6 



L 271 
T 



where 6ps = Ps — Ps(Af)- Eq- ([4.1 6|) is useful to bosonize the Hubbard model in the next 
section. 

V. BOSONIZATION OF THE HUBBARD MODEL 

Bosonization of the Hubbard model can generally be achieved by Tomonaga-Luttinger's 
bosonization techniques ||3^ . However, the existence of spin degree of freedom brings scatter- 



ing processes other than the forward scattering, say the backward and umklapp scatterings. 
They can not be exactly diagonalized. Fortunately in the Hubbard model case, Luther and 
Emery showed that the Hamiltonian still can be diagonalized at particular values of coupling 
constants |^ and the backward and umklapp scatterings can develop gaps for spin- and 



charge- waves respectively. Then reasoning with renormalization group analysis [^, one 



expects that the backward and umklapp scatterings are irrelevant. In this section, we will 
confirm these ideas in the bosonization approach based on Bethe ansatz equations, which 
we have developed in Sec. II to bosonize the C-S model. 



23 



A. Bosonized Hamiltonian and Lagrangian 



Because of the similarity between ( p^.8| ) and ( [4 .21 ), the grand partition function of the 
Hubbard model can be put into a form like (|2.21|). And again this grand partition function 



corresponds to an ideal fermion system with a complicated, T-dependent dressed energy 
ec{k, T). To exploit this fermion representation, in the case of the C-S model we noticed that 
e{k, T) = e{k, 0) + 0(e~l''l/^), so that we could use e(/c, 0) to define an effective Hamiltonian 
to derive the low-T grand partition function. However, this trick can not repeated for 
the Hubbard model, since Scik^T) = Scik) + 0{T'^/u). So we can not simply ignore the 
T-dependent part of ec{k,T) in discussing the low-T thermodynamics of the system. 

Instead of the original thermodynamic potential ( |4.2| ), we begin with the low-T thermo- 
dynamic potential ( |4.16| ). Now, both edk) and es{Sps) are T-independent. In principle, one 
may try to introduce two kinds of fermions with dispersions edk) and es{Sps) respectively. 
However, the fermions with \ps\ > p^p + 6 and \ps\ < PsP ~ ^ are not defined. (One may 
take these states to have zero energy, but it leads to degeneracy deep inside the Fermi sea.) 
Rather, noting the cut-off (5-independence of the low-T thermodynamic potential, we first 
keep a finite 5, and take the limit after all calculations. Thus we rewrite ( [4.16|) as 



fi(T) 




L 


5-0 L 


n{T, 6) 




L 


L L 



(5.1) 



with 



^MH = r dkln{l + e-^^C^'i)/^), 

L 27r J-n 



Os(T, 5) ^ _T_ 
L 271 

Here es{ps, 6) is defined by 



dp,\n{l + e-'^^P^'^^/^). (5.2) 

J —riiTT 



±Vs{PsTPsf), PsF ~ S < \Ps\ < PsF + S, 

£s{Ps, 5) = I ±Vs{k T kF)e~^P^^P^^\/^ , 0<\Ps\<PsF-S (5-3) 

or \ps\ > PsF + ^■ 
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If we keep 6 until all calculations are finished, the degeneracy of the states far from the Fermi 
surface will be removed. Now both edk) and es{ps, S) become T-independent. Therefore, the 
low-T grand partition function can be obtained by an effective Hamiltonian that incorporates 
two kinds of fermions with dispersions edk) and £s(ps, 5): 

k Ps 

= Tre-^-^«ff. (5.4) 
The effective Hamiltonian Hes is given by 

-f^eflf = Y.ec{k)clck + ^es{ps: 5)sl^Sp^ + 5Eo, (5.5) 

k Ps 

where c|. and are spinless fermion creation operators, which can be identified as the 
charge- and spin-excitations. Possible zero-mode excitations are included in the last term, 
SEq, of the effective Hamiltonian. 

Now we are at a position to bosonize the theory, as we have done for the C-S model, 
except that there are two kinds of fermions. Because of the continuity of e'^{k) and e'g{ps, S) 
at Fermi surface, the bosonization is standard. At low-T, only low energy excitations near 
the Fermi surface are relevant. So, instead of the full expressions of the dispersions, the 
linear dispersions are employed 

ef{k) ^ ±e'^{kF){k ^ kp) = ±Vc{pTPf), 

^fiPs) = ^VsiPsTPsp), (5.6) 
and the effective Hamiltonian is separated into the right- and left- moving parts 
-f^eff = Hc,+ + He - + Hs,+ + Hg- + SEq, 

H.,±^T.^fikHc„ Hs,±^J2sf(ps)sls,^. (5.7) 

P Ps 

By using the well-known bosonization technique, we have 



J+) - ^ . J „ . J-) _ V ■ r ■ 

h'c,q ~ • '•-p+q'-^p ■■> Hcq ~ • ^p-q^P ■) 

Ps^q ~ ^ ■ ^Ps+q^Ps PsJ ~ ^ • ^Ps-q^Ps • • (^-S) 

Ps^PsF Ps^-PsF 
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They satisfy the following commutation relations: 



L 
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(±) 1 _ , , 



r (i) (±)t ^ _ ^ . 
[Pc{s),q^ Pc{s),q'i " ^^"9.9'' 



[^cW,±,P;(4J = ±^^c(s)P,,(i),g, (5.9) 

which describe two 1-d free phonons with the sound velocities Vc and Vg- At /i = 0, the 
bosonized Hamiltonian satisfying ( ^.91) is given by 



1 TT 

Hb = Vc ^{alttg + a\aq) + -y{vc,nMI + Vc,.jJI) 

q>0 

q>0 '^^ 



(5.10) 



where 



hg = V2vr/gLp(+), bg = p7i/qLpi~^;i^ (5.11) 

which are normalized boson annihilation operators. The cut-off (5-dependence only con- 
tributes a constant to ( ^.10|) , which vanishes as 5 0. So, we can take 5 — > at this stage, 



which does not affect any results that we will obtain below. 

The bosonized total momentum operator may also be obtained by bosonizing the fermion- 
ized P = Ep cj^p + ^Is^Ps- It reads 



P = q{alaq - aja,) + nJcijicQ + M^/L) 

q>0 

+ E - + ^JsMs/L. (5.12) 

q>0 



In the coordinate-space formulation, the charge- density field and spin-density field are 
j;iven by 

Pc{s){x) = Pc{s),r{^) + Pc{s)A^)^ (5-13) 
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with 



Here Mc{s) = Mc(s),r + Mc(s),l- The conjugate field </)c(s)(a;) of Pc{s){x), which obeys 
[0c(s)(x),pc(s)(x')] =i6{x- x'), is 0c(s)(a;) = (Pc(s),r{^) + <Pc{s),l{^) with 



g>0 V '^(IJ^ 



^ q>0 V '^Q^ 



(5.15) 



where Jc(s) = Jc(s),r + Jc{s),l- Here Mc(s) and Jc(<() are operators with integer eigenvalues 
obeying the selection rule (|3.24| ), and 0co = <Pco,r + 0co,l and 0^0 = (pso,R + 0so,l are angular 
variables conjugated to and M^: [00(5)05^0(5)] = ^■ 

In the coordinate-space formulation, the bosonized Hamiltonian ( |5.1(J| ) reads 



Hb = Hb,c + Hb,s, 

Hb,c = ^ rdx[U,{x)^ + {d^X,{x))% 
Zir Jo 

Hb,s = ^ f'^dxin.ixY + {d,X,ix))% (5.16) 
Ztc Jo 

where Hc(s)(x) = tt^IIs)Pc{s){x) and Xc(s) = ^c{s)'^'^c{s){x). With the Heisenberg operators 
A{x,t) = e^^'^^ A{x)e~^^'^^ , the Lagrangian density is given by 
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C=C, + Cs, (5.17) 

Cc = ^9«Xe(x,t)(9°X,(x,t), (5.18) 
/vr 

Cs = ^d^Xs{x,t)d''Xs{x,t). (5.19) 

ZTT 

The Lagrangians ( p.l8| ) and ( |5.19| ) give rise to two c = 1 Virasoro algebras. We see that there 
is no mixture between spin and charge sectors in the total Lagrangian ( ^.17[ ). Therefore, 



in charge-spin basis, the low energy behavior of the system is characterized by the direct 
product of two independent Virasoro algebras both with central charges c = 1 W^ - 

Note the selection rule (|3.24|) ; then one can take -Dc(s) = 2Jc(s) = Mc(s)mod(l). The 
zero-mode part in Hb,c{s) becomes 

So each sector has the duality relation as in the usual c = 1 CFT. 

Changing the right-and-left-moving representation to the 6-(f) representation with 6 = 
4>R — 4>L, which Haldane has used in discussing the Luttinger liquid theory we see that 
our bosonization theory precisely agrees with the Luttinger liquid theory for the Hubbard 
model proposed by Ren and Anderson The selection rules chosen in [|r^ were 



Ju = Mtamod(2), (5.21) 

where 

Jc=^^, Js = ^^. (5.22) 



The selection rules ( ^.21|) of J| and are a restriction on the quantum numbers more 



constrained than that given by ( |3.24|) . For example, the choice Mc = Mg = 1 ,Jc = 1/2 and 



Jg = —1/2, which equals to M^- = 1, M| = and J| = and Jj^ = 1, is forbidden by ( |5.21| ) 



but is allowed by ( p.24| ). This choice corresponds to a charge-1, spin-up boson excitation 
near kp as we will see below. We believe that ( ^.24| ) gives the correct selection rules for the 
quantum numbers, since it is derived on the basis of the Bethe ansatz equations. 
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B. Correlation functions of single quasiparticles 



In order to calculate the correlation functions in the low energy limit, we first need to 
determine the allowed operator in the theory. Since there are two angular variables, 0co and 
(pso, there are hidden symmetries in the theory under 



'c(s) <Pc{s) 



2n. (5.23) 



— 1/2 

Then Xc(s) is compactified on a circle with radius Rc{s) = X^is-) ' 



Xc(s) ~ Xc(s) + 27ri?c(s)- (5-24) 



Operators V[X] are allowed if only if they are invariant under ( p.23| ) (or ( |5.24| )). Because 



the theory is described by the direct product of two Virasoro algebras, the primary fields 
of the theory have a similar structure. The primary fields satisfy the periodic boundary 
conditions are 



gi{McAy'-Jc/Ay^)Xi,,(x) 



^i(MsXl^''+Js/xV^)XR^,{x) 



^i(Ms\V^-.Js/\V'')XLA^) ■ . (5.25) 
The fermion operator, say of charge-1 and spin-up, is given by = Mg = 1: 



Jc+Js=odd 



where Oj are constant value operators satisfying [Oj, Oji] = i[{Jc + Js)/'^ — {J'c + J's)/'^]- 



The above fermion operator satisfies our selection rule ( |3.24D and, since Jc + Js = odd, 
the selection rule ( ^.21| ) of Ren and Anderson as well. However, our rule (|3.24| ) allows 
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more possible quantum numbers, say charge-1, spin- up bosonic excitations described by the 
primary fields: 

Jc+^s=even 

-Js/\V'^)Xl^^{x) . ^ (5.27) 

The selection rule for bosons reads 

JTa = (Mt,; + l)mod(2), (5.28) 
which is different with that for fermions, (|5.21|) . 



Besides the fermion and boson excitations, there are two kinds of anyonic excitations 
which are charge-1 but non-periodic primary fields. The anyon operators read 

iT/t _ iTrt p«Ay^(XB,c-^i,c) 

= ^t^e*^^^'(^«'^-^^'^\ (5.29) 

which obey 

^1 (x)^\ (x') - e^^-w^^l (x')^l (x) = 0. (5.30) 

Thus, the statistics of is interaction-dependent, while \[^t^^ is a semion, in the half way 
between fermion and boson. 

The dynamic correlation functions of fermionic and bosonic single-particle operators, 
( p.26| ) and ( |5.27| ), can be easily calculated: 

G,.(x,t) = (*^(x,t)^U0,0)) 



oo 

ra=— oo,m=4n+l 



2 trriKp 



oo 

+ E 

?i=— oo,m=4n+3 

a^R,c ^l!c x^je*, (5.31) 
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GB{x,t) = (^B(x,t)^MO,0)) 

oo 

^ Cm 

ra=— oo,m=4ra+l 

oo 

ri=— oo,m=4n+3 

where A^, B^, Cm and Dm are regularization-dependent constants, and x/j_c(s) = x — fc(s)^ 
and xl,c{s) = x + Vc{s)t. e"^'^ = 2Ac. To derive (|5.31| ) and (|5.32|) , we have taken the leading 
term only for each given m. For the fermion case, is taken to be ^ or —| with respect 
to m = 4n + 1 or 4n + 3. For the boson case, Js is taken to be —| or | with respect to 
m = 4n + 1 or 4n + 3. We see that the only difference between fermion and boson is that 
when the fermion has a left(right)-moving spinon-part in a given mfciT'-oscillation, the boson 
has a right (left)-moving one. They have the same momentum distribution singularities: 

n{k) ~ const. - const. sgn(A; - kF)\k - kpl^^''^^''^-^^^^ (5.33) 

for /ciT'-oscillations and 

n{k) ~ const. - const. sgn(fc - 3kF)\k - 3kF\''^''~'^+''^-^'>/\ (5.34) 



for 3fci;'-oscillations, etc. They are consistent with those given in In the strong coupling 
limit, U/t oo, one has (e'^^ + e^^ - 2)/4 = | and (Oe'^^ + e^^ - 2)/4 = | as e"^^ = |. 
This implies that there is a similarity of boson and fermion in the two-component Luttinger 
liquid theory. 

The anyon correlation functions read 

CxAx,t)^{^xA^,t)^iMO)) 



n= — 00,771=471 + 1 

- i (me*' )2 - i imkp 
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n=— oo,m,=4n+3 



(5.35) 



and 



GA.(x,t)^(vI/,^(x,t)vI/t^(0,0)) 

^m-^R,c -^Lfi •^R,s-^L,s^ 

n=— oo,m=4n+l 



8 irriKp 



oo 



+ E B>S""'"""''xZf--""''x3xZJe-^. (5.36) 

n=— oo,m=4n+3 

The momentum distributions corresponding to (|5.35|) vanish faster than those for 
fermions or bosons near mkp. So, this kind of excitations may be more difficult to be 
observed. The distributions corresponding to ( p.36| ) near (4n + l)A;ir, say kp, also decay fast 
but those near {An + 3)A;f, say 3fci?, do decay slower: 

n{k) ~ const. - const. sgn(fc - 3fcir)|fc - 2>kF\^'''~"'^''^-^^'^ 

7 

^ const. — const. sgn(fc — 3A;ir)|fc — 3fci?| 8, U/t ^oo. (5.37) 

So, the semion may possibly be observed. 

As shown by ( p.46| ), for the C-S model, it is the same density operator no matter that 
the single-particle excitations are fermionic, bosonic or anyonic. This remains true for the 
charge density and spin density in the Hubbard model. And so for the charge-charge and 
spin-spin density correlation functions. We can define the pairing operators in the theory. 
One finds that the most important contributions to the spin-singlet and triplet pairings, 
up to 2A;i?-oscillation, still come from the fermionic (or bosonic) excitations. All of these 



correlation functions have been calculated in the reference , and can be reproduced easily 
within our approach. 
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VI. CONCLUSIONS 



We have developed a new approach for bosonizing 1-d exactly soluble many-body mod- 
els at low energy. In this approach, one first derives the low temperature partition function 
from the Bethe ansatz equations, and then bosonizes the partition function. The asymp- 
totic correlation functions can then be calculated by field theory techniques. Because our 
bosonization started with the Bethe ansatz equations, we call it the bosonization based on 
Bethe ansatz equations. We explicitly carried out this program with the C-S model and the 
Hubbard model as examples. Some new results were obtained during this exercise. 

For the C-S model, we showed that the low energy effective theory is described by 
two different classes of c = 1 CFT's, depending on whether the Hamiltonian describes the 
interactions between bosons or between fermions. Using the standard terminology in CFT 
3^, it is a c = 1 CFT compactified on a circle 5*^ for the bosonic case, compactified on an 



interval jZ^ for the fermionic case. A piece of evidence for this difference is that the charge- 
1 primary field in the two theories are different. It is a fermion in the theory of interacting 
fermions, while it is a boson for the interacting bosons. Their correlation functions have 
totally different asymptotic behavior. Also our bosonization of the C-S model prepares some 
techniques for bosonizing mult i- component models. 

The repulsive Hubbard model at less than half filling is an example of the two-component 
models. In the above the spin-charge separation in this model has been explicitly verified 
not only for the strong-coupling limit but also for finite U. The low energy theory is shown 
to be controlled by the direct product of two usual c = 1 Virasoro algebras. We confirmed 
this point by the finite-size scaling in the temperature direction. Furthermore, the explicit 
low-T grand partition function enabled us to identify the system with a two-component Lut- 
tinger liquid theory by using the bosonization technique. We showed that the backforward 
scattering processes are suppressed at low energy, in agreement with the renormalization 
group analysis. We pointed out that not only the fermionic excitations are allowed near 
the Fermi surface, but also the bosonic and anyonic ones. The single-particle correlation 
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functions of those excitations were systematically calculated in a simple way. 

Our bosonization approach can be easily generalized to any 1-d many-body models that 
is exactly soluble by Bethe ansatz, including thermodynamic Bethe ansatz, and provides a 
simple way to calculate asymptotic correlation functions using field theory techniques. 

This work was supported in part by the US NSF grants PHY-9309458 PHY-9970701 and 
NSF of China. 



APPENDIX A: FINITE SIZE SCALING IN THE L-DIRECTION FOR THE C-S 

MODEL 

This appendix is devoted to resolve a puzzle raised in where it was found that if 
we take k = ^j^n in the thermodynamic limit, naively the discrete sum for the total energy 
of the C-S model differs from the continuous integration by — This seems to imply a 



central charge c = A, in conflict to ( p.l9|) . However, as we see below, a careful consideration 



for the finite size scaling shows that the naive result is wrong. There is an additional term 
of order 0(1/ L^) when we go to the continuous limit of ( |2.5| ) from the discrete version. This 
will cause a rescaling of the particle density together with that of the total energy. Taking 
into account this subtlety, we will finally achieve the correct finite size scaling of the ground 
state energy 

Eo,L-E, = -^ + Oi^). (Al) 
To see this, a useful formula relating the discrete sum to the integration is 

1 ^2 T r{N2-l/2)/L 

L nt^^ L 7(7Vi+l/2)/L 

+ ^[/m-i/2)/L) 

- n{N, + l/2)/L)] + 0{l/L'). (A2) 



Using (|A2|) , the discrete version of the density piik) can be written as 
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PL{k) = - + {l-\)pL{k) 



+ ^(1 - \)\^[5{k - kp) - 5{k + kp)]. 
akp 



12L2 



(A3) 



The discrete version differs from its continuous counterpart by a term of the order 0{1/L'^). 
Denote 



Piik) = p{k) +pi(/c); 



(A4) 



one has 



p.W = -^(i-A)[J;*(* + M-5^*(*-M]. 



(A5) 



while p{k) = 1/2'kX. Now let's examine the finite size scaling of the ground state energy. 
Using ( |A2|) again, we have 

— = / dkpL{k)eo 
L J-kp 



dkp{k)to{k) + dkpi{k)eo{k) - ——. 

-kp J—kp Uij 

The second term of the last equation is easy to calculate and one has 



'^^dkp,{k)eoik) = ~^{l-\). 



Substituting this into ([A6D , we have 

EoL Eq hvf 



6L2 



This confirms the result of (|2.2CI|) . 



(A6) 



(AT) 



(A8) 



APPENDIX B: PROOF OF e{T) = e(0) + 0{T'^/v) 



In this appendix, we derive the low temperature expansions of Scik, T) and ^^(A, T), by 
solving the integral equations (^4.3|) by iteration. 



35 



At the zeroth-order we take 

e^^Hk,T)=e^^\k), .f(A,T) = .r(A) = ^. (Bl) 
Then at the first order, i.e. after one iteration, the dressed energies are 

£«(A;,T) = 40) + O(Te-'^/2T)^ (B2) 

eW(A,T) = 4°)(A) + - y_J„^ dkK,{X - sinfc)4°)(A;) 

-- r rfA;Ki(A-sinA;)ln(l + e-"'°^) 

+0{Te-^'^^), (B3) 

where k^p^ is determined by s'f\k) = 0. We denote s'^^^k) = e''^\k) and the first three terms 
on the right hand side of ( P^ ) by £:^^^(A). The fourth term on the right hand side of ( [B3| ) 
is proportional to T'^/Vc, if we consider only the contributions from near the Fermi points. 
Thus, we have 

4^)(A;,T)=.W(fc) + 0(Te-'^/^^), 

4'nA,T)=e«(A) + 0(TV^;,). (B4) 

To get 4^)(A;,T) and ef\k,T), it is enough to replace e^P{k,T) and e[^\k,T) by 4^)(/c,0) 
and s'^p{k^ 0) on the iteration equations. Then by using similar techniques, we have 

ef(A;,T) = ef(A;) + 0(rVt;,), 

ef\\T) = s^P{X) + 0{T'/u). (B5) 

where s^^\k) and £^^^(A) are given by the iteration of integral equations ( |3.7|) and v could 
be Vg or Vc- Repeating the iteration, we finally have 

lim ei'^Hk^T) = lim ^"^(A;) + OfrVz/), 

Jim 4")(A, A;) =jim4")(A) + O(TVz^). (B6) 

36 



Because of the convergence of the series for e^^l~^{k,T), ( p6|) imphes that 



e,{k,T) = e,{k) + 0{Tyu), 
es{X,k) = es{X) + 0{Tyu), 



(B7) 



which is what we intend to prove. 



APPENDIX C: AN ALTERNATIVE PROOF OF ( gTT2D 



We give an alternative proof of (|4.12| ). For the simphcity, we consider the h = case in 
which, as Lieb and Wu pointed out, there are no particle-hke spinon excitations. The total 
energy is given by 

EiT) 



dkp,{k,T)ef\k,T) 

/oo 
ciAp,(A,T)4°)(A,T). 
-oo 



Substituting 



3[) into ( |C1|) , one has 

E{T) 



+ 



+ 



dkpc{k,T)ec{k,T) 



-ea{k,T)/T- 



cos k 

271 



r dkT\n[l + e 

J —TT 

/OO 
dXp,{\,T)Ki{\~smk) 
-oo 

/oo 
dXps{X,T)e.,{X,T) 
-OO 

/oo 
rfArin[l + e-^=(^'^)/^] 
-oo 

X [— dkpc{k,T)Ki{X-smk) 
.27r J-TT 

dX'psiX',T)K2iX- X' 



1 r°° 
2n J-oo 



(CI) 



(C2) 



Then divide the integration range of k into the intervals [— vr, —kp], [—kp, —kp + 5], [—kp + 
(5, fci? — 5], [fci? — 5, fci?] and [fci?, vr] and the interval of A into [—Xp, —Xp + 5], [—Xp + 5,Xp — 5] 
and [Xp - 5,Xp\. Use the fact that e'l^^'^'^)'/^ decay rapidly as T ^ and (|3ll|) . After 
some algebras, one gets 
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E(T) 1 ri^F 



27r J-k 



F 



\^c{k)\ 



rkp+S tuwirr 1 



+2 / " dkT\n{l + e-l^=('=)l/^)(p,(fcp) - — ) 

\es{X)\ 



f-A 



+2/^ dATln(l + e-l^=WI/^)p,(Ai.). (C3) 
Here 5 is a cut-off. Using the integral formulas 



dx- 



h 1 + 
we have, in the limit T ^ 0, 



|^°°rfxln(l + e--) = ^, (C4) 



E(T) 1 1 r'^p 

= / rfA;£,(A;) + — / dki,{k, T) 

7rT2 7rT2 vtTV. 1 



On the other hand, ai h = 0, in terms of ( [4.4D , we have an integral equation for ec{k,T), 

e^{k,T) = -^f{k), (C6) 

/(A;) = cos A;^K(sin k — sin fc^) 
1 ff^F 



+— f " dk' cos k'K {sin k- sin k')f{k' 
2tt J-kp 

/kp dj^ ^ 
— ec{k,T), one has 
-kp 271 



F 

dk 

kp 

f^F dk ^, ttTV 1 \ ^^^^ 



So 



EiT) 1 fkF vrT^ vrT^ 

- ' c/te,(A;) + ^ + ^ (C8) 



L 271 J-kp 6fs 
Finally, using the thermodynamic relation 



EiT) ^ m. (C9) 



we obtain the thermodynamic potential as given by (f4. 12|) 
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